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Abstract 

A Green's function formalism is used to calculate the spectrum of localized modes 
of an impurity layer implanted within a ferromagnetic thin film. The equations 
of motion for the Green's functions are determined in the framework of the Ising 
model in a transverse field. We show that depending on the thickness, exchange and 
effective field parameters, there is a "crossover" effect between the surface modes 
and impurity localized modes. For thicker films the results show that the degeneracy 
of the surface modes can be lifted by the presence of an impurity layer. 
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1 Introduction 



The spectrum of elementary excitations in solids is known to be modified 
by the presence of surfaces in the medium in which they propagate. This 
occurs because an interface breaks the translational symmetry of the medium 
and can also modify the microscopic scheme of interactions due to surface 
reconstruction effects. The modifications of the spectrum can be signaled by 
the appearance of surface localized modes, besides the volume excitations. 
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The translational symmetry of a crystal can also be broken by the presence 
of defects or impurities in the lattice, and localized excitations may result. 



Localized excitations of magnetic media have been widely studied, theoret- 
ically as well as experimentally [1] and several models have been proposed 
to elucidate the dynamics of surface modes in a variety of magnetic systems. 
Among these models, the transverse Ising model has been shown to give a 
good theoretical description of real materials with anisotropic exchange (e.g., 
C0CS3CI5 and DyPO/i), of materials in which the crystal field ground state is 
a singlet [2], and has also been employed as a pseudo-spin model for hydrogen- 
bonded ferroelectric materials [3]. Previous calculations have obtained the spin 
wave (SW) spectrum of Ising ferromagnets in a transverse external field for 
semi-infinite systems [4] and films [5]. These results were obtained for ideal 
media, i.e. structures without defects or impurities. Studies [6,7,8,9,10] con- 
cerning excitations of semi-infinite systems described by the transverse Ising 
model have predicted the existence of localized SW modes, as well as bulk 



Within the quasi-particle approximation, the effect of a localized perturbation 
in a periodic structure can be worked out exactly using Green's functions tech- 
niques. These techniques have been extensively used to calculate the impurity 
modes in infinite and semi-infinite Heisenberg ferromagnets and antiferromag- 
nets [11,12]. 

The purpose of this paper is to present the first calculation of the SW spectrum 
of a thin ferromagnetic film with an impurity layer in the framework of an 
Ising model in a transverse field. Dispersion relations for the bulk, surface and 
impurity SW modes are obtained by finding the poles of the Green's functions 
associated with the spin operators of each layer in the film. This calculation 
generalizes a previous study of the spin dynamics of an impurity layer in a 
semi-infinite ferromagnet, presented in Ref. [10]. 



2 Green's function formalism 



Let us consider a retarded Green's function, defined in terms of two time- 
dependent operators A(t) and B(t), in the Heisenberg picture, as: 



SW. 
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which, in frequency representation, can be shown to satisfy the equation of 
motion [13], 



u{{A;B)) u = ±{[A,B]) + {{[A,7{\;B)) u , (2) 



where a; is a frequency label, and Tt is the Hamiltonian of the system. In the 
present case we employ the Ising Hamiltonian in a transverse field 

n = -\T. JaS^-h^Sf, (3) 

<i,j> i 



where Sf{a = x,y,z) is a spin operator component at site i, h is the applied 
transverse magnetic field, taken as hs in the surface and hj in the impurity 
layer, and Jij is the exchange coupling between nearest-neighbor sites % and j. 
The exchange parameter is taken as J in the bulk region, as Js in the surface 
and as Jj in the impurity layer. The exchange coupling between a spin in the 
surface layer or in the impurity layer and a spin in the adjacent bulk layers is 
assumed to be J' s and Jj, respectively. In the special case where the impurity 
layer and the surface are adjacent we use J$i as the exchange interaction. 
Next, we construct the equation-of- motion 



+ 

[S?,H];Sl 



(4) 



where / and m are labels assigned to each lattice site. Thus we obtain an 
infinite chain of coupled equations, which can then be decoupled using the 
Random Phase Approximation [14], 



[S?S];S&)) i n(£?)((S];S5 i )) u + (5) 



where I ^ j. Using this approximation and inserting the Hamiltonian Eq. (3) 
into the equation of motion Eq. (5), we obtain 



A(u>)((Sr,SM) u =hK! 5 £- (6) 
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where 
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A{u) 



- E JuRj ( E W) - h2 

3 V i / 



with Rf = (Sf) and i2f = (Sf). 

Equation (6) may be solved to obtain the spin-dependent Green's functions at 
any temperature. In this work we focus on the paramagnetic phase, at T > T c , 
where T c is the Curie temperature of the ferromagnet. In this regime, due to 
the absence of net magnetization, we have (<Sj\ = 0. For spin S=l/2, we 
use the static spin average in the x direction calculated from the mean-field 
theory [3] 



(7) 



where, depending on the site layer , the index £ stands for S, B, I, for surface, 
bulk, and impurity layers, respectively and ks is Boltzmann's constant. After 
some straightforward algebra, the Green's function is found as: 



G 



nnl 



2ir 



hm 



uj 2 - h 2 +4/ii?f J 7 (q||) 



(8) 



where n and n' are layer indices and 7(qn) is the structure factor, given by, 



t(Q||) = ^[cos(g x a) +cos(g w a)]. 



(9) 



Some mathematical techniques, including general recursive algorithms for lay- 
ered systems, are available to solve these equations. Here we obtain explicit 
results using an approach analogous to earlier calculations for pure ferromag- 
nets [8]. We write the coupled equations in a matrix form 

(A + A)G = -(-L)l (10) 



where A is a iV x iV tridiagonal matrix, with N being the total number of 
layers in the film. The elements of A are defined by 

[A)]z,m = d 8 hm - (5/_i, m + 5 l+ljTn ), (11) 



which represents the system without impurities, and A is a iV x N matrix con- 
taining the information regarding the perturbing effects of both the surfaces 
and the impurities, with elements given by 
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Ai,i = A N<N = A s (12) 
^i,2 = A 2j i = Atv,at-i = Ajv-i,jv — 
A M = Aj 

Aj+i,/ = Aj_ 1; ; = A i)J+1 = A^_! = A IB , 

where the index / refers to the impurity layer, and the other elements are set 
to zero. The parameter d is defined by: 

^-fr 2 + 4/^J 7 (q||) (u) 
hRU ' 1 ' 



The parameters A 5 , As#, A/ and A/ B , which depend on the surface and on 
the impurity are 

= y-^^rA^ _ d (14) 



k5B = — — 



(15) 



y-^^q,)., (16) 



and, 



A„ = (17) 



Finally, we rewrite the Green's function in matrix form: 
1 



G = - 



I + Ao^A)" 1 Ao 1 . (18) 



,2ttJ, 

The inverse of the tridiagonal matrix A can be found in the literature. 
3 Surface and impurity modes 



In order to calculate the dispersion relations of the localized SW modes, one 
has to find the poles of the Green's function Eq.(18). These are obtained by 
solving the determinantal equation [8] 

det(I + A _1 A) = 0. (19) 
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The elements of the inverse matrix Ao 1 involve a factor x, which is defined 
in terms of the diagonal element of A as 

d = x + x~ 1 . (20) 



Solutions of Eq. (19) involving \x\ = 1 give us the frequencies of the quantized 
bulk modes. Such solutions correspond to x — exp(iaa), where a is a set of 
real discrete values (0 < a < n) which depends on the film thickness and 
can be obtained from the determinantal equation. Therefore, from Eq.(13) we 
obtain the relationship: 



u B (q\\, a) = {h 2 - 2hR x Jry[cos(q x a) + (21) 

cos(q y a) + cos(o;a)]}2 

where qy is a two-dimensional wave vector. 

We are interested in understanding how the surface and localized impurity 
modes affect each other. In order to do that, we start our analysis by consid- 
ering a film with 5 layers, where the first layer is taken as the impure one. In 
this case, the A matrix can be expressed as: 



A 7 


A/b 








































A 5 








OAs 


Ai 



(22) 



Figure 1 shows a plot of frequency as a function of wave vector q x a/n, where 
q y = 0, for an impurity layer located on the surface of a 5 layers film, with 
the parameters: R x = 0.25, Jj = 1.5J, j\ = 1.25J, J s = 1.5J, h s = 1.5/i, 
hi = 1.35/i and J SI = 1.0J. The shaded area in the graph corresponds to the 
region containing the volume modes of the infinite case. One can observe the 
presence of two surface branches located above the bulk region. The branches 
are associated with the two surfaces of the film, one of these being the impurity 
layer. The figure also shows a "crossover" effect (with mode repulsion) between 
the surface and the surface impurity modes at q x ~ 0.22. The existence of this 
crossover is dependent on the strength of the impurity exchange, as can be 
seen in the next graph. The results in Fig. 2 were calculated for a smaller value 
for the exchange constant in the impurity layer, namely Ji = 0.5J, with the 
remaining parameters being the same as those in Fig. 1. The graph shows two 
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Fig. 1. SW dispersion for a thin ferromagnetic film with an impurity layer at the sur- 
face, for J j = 1.5 J. The bulk band is represented by the shaded region. Parameters 
in the text. 



localized SW branches. The lower frequency mode does not display any shift 
in relation to the result in Fig. 1. On the other hand, the upper SW branch 
shows a significant shift in comparison with the corresponding branch in Fig. 
1. Thus, we may associate the lower frequency branch with the pure surface of 
the film, whereas the upper branch can be associated with the surface impurity 
layer mode, since the shift occurs as a direct consequence of the smaller value 
of J/. Also, due to the frequency shift, no mode repulsion effects are found. 



Let us now consider the situation in which the impurity layer is located in the 
interior of the film. Results for this case are shown in Fig. 3, for a 5 layer film, 
in which the impurity layer corresponds to the 3rd layer of the film. The A 
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matrix for this configuration is: 
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(23) 



The SW frequencies are obtained following a procedure similar to the previous 
case. The parameters are hi = lA5h, Jj = 1.5J, Jj = 1.25J, hg = 1.35/i, Js = 
0.5J, J s = 1.0J, Jsi = 1.0J, with R x = 0.25. The plot shows three coupled 
modes. For small wave vectors, there are two degenerate frequency branches 
related to the surfaces, as well as a lower frequency one (impurity mode). 
As the wave vector increases, the degeneracy is lifted. In the range 0.57T < 
q x a < 0.77T, there are three distinct modes. That happens as a consequence of 
a double crossover, due to the mixing of the surface and impurity modes. For 
larger values of q x , the surface modes are again degenerate. 

Figure 4 shows a graph of the frequencies of localized SW modes as a function 
of the effective field hi in the impurity layer. We used a film with the same 
parameters as the one discussed in Fig. 3. Results are shown for three different 
values of in-plane wave vector, namely q x a = O.Oln (solid line), q x a = 0.57T 
(long-dashed line) and q x a = l.On (dotted line). In contrast with the semi- 
infinite case presented in Ref. [10], the graph displays a rich spectrum of 
localized modes. The discontinuity in the graph is a consequence of the merging 
of some of the modes with the bulk band for certain values of hj. For each 
value of wave vector, the graph shows two different behaviors of the localized 
SW branches: some branches display a strong dependence on hi, whereas the 
remaining branches are not noticeably influenced by a varying field, except at 
larger values of hi, where mode mixing effects are observed. These contrasting 
behaviors reflect the distinct nature of these modes, with the impurity modes 
being associated with the branch that is strongly influenced by the hi field. 
These impurity modes have acoustic nature for small values of hi. As the 
effective field is increased, the impurity mode merges with the bulk band 
and eventually reappears as an optical mode. For the parameters used in this 
calculation, the surface modes are found above the bulk band. 
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Fig. 2. SW dispersion, for a system similar to the film of Fig. 1, with Jj = 0.5 J. 



This difference in behavior is also evident in the results of Fig. 5, which shows 
localized SW frequencies as a function of the effective surface field hs, for 
the same wave vector values as in the previous results and hj/h = 1.45. The 
graph shows the influence of hs on the surface modes. For small values of this 
effective field, the surface modes are degenerate and have acoustic behavior. 
On the other hand, the impurity modes are not affected by the variation of hs- 
As the field is increased, the surface modes merge with the bulk band, and an 
acoustic impurity mode is observed. Mode mixing effects are also found. For 
larger values of hs, surface modes are found above the bulk band and there 
are two impurity modes: one optical branch, and one acoustical. 



4 Conclusions 

In summary, we developed a Green's Function formalism for calculating the 
SW spectrum of a ferromagnetic thin film with an impurity layer, in the 
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Fig. 3. SW dispersion, for an impurity layer in the interior of the film, correspond- 
ing to the third atomic layer. The bulk band is represented by the shaded region. 
Parameters in the text. 

framework of the transverse Ising Model, using the RPA to obtain equations 
of motion for the bulk, surface and impurity SW modes. We calculated the 
Green's functions analytically through an inhomogeneous matrix equation in 
a closed form. The formalism applied for any temperature, and we 

obtained explicit solutions for T > T c . The results showed the effect of the 
exchange coupling between spin sites in the impurity layer on the localized 
modes of a film in which the impurity layers correspond to one of the surfaces. 
We also obtained results for the effect of the position of the impurity layer in 
the film. We found that the placement of the impure layer has a strong effect 
on the spectrum of localized modes. In all cases, we observed strong mode 
mixing effects between localized impurity and surface modes. In addition, we 
calculated the localized SW frequencies as functions of the effective fields hi 
and hs- This allowed us to analyze the nature of the modes, by observing 
the response of the frequency branches to the variation of the effective fields. 
In the limit N — > oo, the results agree with previous calculations for semi- 
infinite systems [10]. This formalism can also be applied to obtain the spectral 
intensities of the bulk and localized SW modes. Currently, we are extending 
the formalism in order to include the effects of localized impurities in thin 
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Fig. 4. Localized SW frequencies as a function of impurity magnetic field for three 
values of wave vector. Parameters in the text. 

films. 
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